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We analyze the viscous damping of r-mode oscillations of compact stars, taking into account non- 
linear viscous effects in the large-amplitude regime. The qualitatively different cases of hadronic 
stars, strange quark stars, and hybrid stars are studied. We calculate the viscous damping times 
of r-modes, obtaining numerical results and also general approximate analytic expressions that ex- 
plicitly exhibit the dependence on the parameters that are relevant for a future spindown evolution 
calculation. The strongly enhanced damping of large amplitude oscillations leads to damping times 
that are considerably lower than those obtained when the amplitude dependence of the viscosity is 
neglected. Consequently, large-amplitude viscous damping competes with the gravitational insta- 
bility at all physical frequencies and could stop the r-mode growth in case this is not done before 
by non-linear hydrodynamic mechanisms. 



I. INTRODUCTION 

A compact star is one of the most stable forms of mat- 
ter in the universe. The only instability that threatens 
its existence is collapse into a black hole, triggered by un- 
stable radial oscillation modes that push the star below 
its Schwarzschild radius. However, there are other unsta- 
ble oscillation modes, the so called r-modes Pll?, which 
damp the rotation of the star by the emission of grav- 
itational radiation [3]. In slowly-rotating stars r-modes 
themselves are damped by bulk and shear viscosity, but 
at high rotation frequencies they are unstable and grow 
exponentially. In a companion paper |1] we analyzed 
these instability regions in detail. The main result of 
that study was that these regions vary greatly between 
qualitatively different classes of stars containing distinct 
phases of strongly interacting matter, but are extremely 
insensitive to unknown quantitative details of the equa- 
tion of state and the transport properties within a given 
class. Therefore, a proper understanding of the r-mode 
dynamics could in the future provide robust signatures 
for the presence of exotic phases in compact stars. 

Since exponentially growing r-modes will destroy the 
star if their growth is not stopped by some non-linear 
mechanism, the fact that fast spinning compact stars are 
observed suggests that such a non-linear damping mecha- 
nism is indeed present. More importantly, even if stopped 
at a finite amplitude, r-modes still strongly emit gravi- 
tational waves and could provide an extremely efficient 
mechanism for the spin-down of compact stars [5] and an 
interesting signal for terrestrial gravitational wave detec- 
tors. Spin-down due to r-modes could explain the ob- 
served absence of fast-spinning young stars despite the 
fact that their creation during a supernova could natu- 
rally lead to a fast spinning remnant. For spin-down via 
r-modes the size of the saturation amplitude is crucial. If 
the amplitude a is too low, it takes too long to spin down 
the star; if the amplitude is too large, a> O (1) , the r- 
mode would disrupt the star's structure, and even before 
this point the r-mode could be destroyed. If the r-mode 
amplitude saturates at an intermediate value, a fast spin- 
down is possible. Previously, various mechanisms for the 



large- amplitude behavior of r-modes have been suggested 
They include the coupling between different modes 
[71 H], the decay into daughter modes and the eventual 
transformation of the r-mode energy into differential ro- 
tation [9,|10J, friction between different layers of the star, 
and surface effects in the star's crust 

Here we study an alternative mechanism that does not 
involve such complicated non-linear dynamical or struc- 
tural effects. It is present already in a standard hy- 
drodynamical description and exploits the fact that at 
large amplitudes the damping due to bulk viscosity in- 
creases dramatically with the amplitude [12H15| . In this 
suprathermal regime, where the deviation from chemical 
equilibrium /i^ fulfills /ia ?J T, the viscous damping could 
overcome the initial gravitational instability and saturate 
the r-mode. However, as shown in [12 , the bulk viscos- 
ity has a maximum as a function of the amplitude and 
decreases again at even larger amplitudes. If the ampli- 
tude exceeds this critical value then the r-mode growth 
cannot be stopped by viscous damping and other non- 
linear dynamic effects [SHTU] are required to saturate it. 
Nevertheless, we find that over a significant region of the 
parameter space the suprathermal enhancement is indeed 
sufficient to saturate the r-mode at a finite amplitude and 
the r-mode can then efficiently spin down the star. This 
is in contrast to certain non-linear hydrodynamical ef- 
fects where the r-mode could completely decay [10 and 
would not be able to cause an appreciable spin-down of 
the star. 

For r-modes with amplitudes sufficiently below the 
maximum, we give general analytic expressions for the 
suprathermal damping time valid for various forms of 
dense matter. For r-modes with arbitrary amplitudes, 
where an analytic evaluation is not possible, we give a 
general expression for the bulk viscosity damping time 
that includes the complete parameter dependence re- 
quired for the analysis of the star's evolution, encoded in 
a two-parameter function that can be numerically com- 
puted and tabulated for different star models. This of- 
fers an explicit framework for the consistent inclusion of 
the r-mode saturation into a star evolution analysis and 
supersedes previously necessary model assumptions 0. 
We will analyze the same star models as in our recent 
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companion paper U and thereby extend this study to 
the suprathermal regime. These include neutron stars, 
hybrid stars and strange stars. In addition, motivated 
by the recent observation of a 2Mq star |16l [T7] we only 
study models that also yield heavy stars. Moreover, simi- 
lar to the analytic results in [4_ , we give approximate ana- 
lytic expressions for the maximum saturation amplitude 
that exhibit the detailed parameter dependence on the 
equation of state and the transport properties of dense 
matter. In addition to the standard fundamental to = 2 
r-mode we also study the saturation of higher multipoles 
and find that they saturate at significantly lower am- 
plitudes. In this work we concentrate on the effects of 
suprathermal bulk viscosity on the damping times |18| 
and the suitably extended concept of the instability re- 
gions, and defer a comprehensive analysis of the star evo- 
lution to future work. 



II. DAMPING OF STAR OSCILLATIONS 

In this section we discuss the prerequisites needed for 
the analysis of compact star oscillations and their damp- 
ing. In particular we discuss the large amplitude en- 
hancement of the bulk viscosity which will provide a 
mechanism for the saturation of r-modes. This topic had 
been studied in [T? and we refer the reader to this work 
for further details. The analysis further requires the sta- 
ble equilibrium configuration of the star and the r-mode 
profile. These have been discussed in more detail in our 
companion article on the damping of small amplitude r- 
modes |1] and we will here only briefly recall these results. 



A. Suprathermal viscosity 

The bulk viscosity describes the local dissipation of en- 
ergy in a fluid element in one cycle of compression and 
rarefaction, driven by some oscillation mode of the star. 
The integration over the whole star yields the correspond- 
ing total energy dissipation of an oscillation mode as will 
be discussed in section |III[ Recently the bulk viscosity of 
large-amplitude oscillations has been studied in detail in 
|12| . which builds on the classic work [14 and provides 
general expressions valid for various forms of matter and 
arbitrary equations of state. 

Bulk viscosity is generally induced by slow weak- 
interaction processes whose rate takes the parametric 
form 




(1) 



where T is the temperature and /^a represents the quan- 
tity that is driven out of equilibrium due to the oscilla- 
tions and its re-equilibration leads to the bulk viscosity. 



The latter is given by the difference of the sums of chemi- 
cal potentials of the particles in the initial and flnal state 
of the relevant weak process. The Xj &re coefficients that 
characterize the non-linear (large amplitude) contribu- 
tion to the re-equilibration rate. The series terminates 
at a finite order N which is determined by the number 
and type of particles in the initial and final states of the 
dominant re-equilibration process and is connected to the 
temperature dependence via 5 = 2N . 

As has been shown in |T2 the bulk viscosity can be 
written in a general form that expresses its full underlying 
parameter dependence in terms of the coefficients d of the 
driving term and / of the feedback term 
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in the differential equation that determines the oscilla- 
tion /iA- These expressions depend on the oscillation 
frequency w, the conserved number density fluctuation 
An/n, as well as the susceptibilities 
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with respect to the density n and the fraction x of a 
particular particle that is driven out of equilibrium. The 
general result can be parameterized in the form [12J 



where the global maximum of the viscosity 
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is taken at 



= ,/ = 1 =^ Tmax = 



TB 



(4) 



and the non-trivial dimensionless integral 
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has to be tabulated for each considered form of matter 
with a particular dominant weak re-equilibration process. 
Since both d and / depend on the temperature and the 
required parameter regions are therefore non-trivial, see 
|19] . this is most conveniently done in terms of the new 

variable d = d f^^^ for the function I {^d, =2{d,f). 

This form of the bulk viscosity is crucial for the study of 
damping times below, which involve an integration over 
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the whole star, and therefore requires the complete den- 
sity and amplitude dependence of the viscosity in addi- 
tion to the temperature and frequency dependence which 
are used in standard analyses. 

In general the bulk viscosity features three distinct 
characteristic regions. Almost all previous analyses of 
r-mode damping have been limited to the suhthermal 
regime /za <C T, where /iA is linear in d ~ An, so the 
viscosity is independent of the amplitude and has the 
analytic form 
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Yet, since the r-mode is unstable and rises exponentially 
it eventually reaches the suprathermal regime /^a ^ T 
which will be studied below. This regime is divided 
into two qualitatively different parts. In an intermediate 
regime /iA is still linear in d but the viscosity strongly 
rises. For / ^ 1 this part allows an analytic solution 
given by 
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which can be combined with the subthermal result to 
give an approximate analytic solution for both regions 
C~ = C<+0(T— T'„iQa;)C~. At even higher amplitudes the 
rise of /iA becomes weaker due to non-linear saturation 
effects so that the viscosity has a maximum and decreases 
again. In the asymptotic limit /ia 3> T , the viscosity 
scales as 



c 



(8) 



As discussed in appendix|B] for the special case of strange 
quark matter an approximate analytic result is possible 
that includes this large amplitude regime. 

In contrast to the bulk viscosity, the shear viscosity of 
dense matter is independent of the frequency and ampli- 
tude of an external oscillation and over certain temper- 
ature ranges its dependence on temperature is approx- 
imately a simple power law. Shear viscosity becomes 
large at low temperatures and therefore it is the domi- 
nant process for damping of the r-modes of cooler stars. 
Correspondingly we can be parameterize the shear vis- 
cosity as 



(9) 



by simply factoring out the temperature dependence with 
exponent a. 



B. Star models 

We study in this work the same model examples of 
compact stars as in our companion paper These in- 
clude neutron stars, strange stars [50^, and hybrid stars, 
and we consider in each star with a standard mass 

of 1.4 Mq and a heavy star with a mass of 2 Mq. For the 
neutron stars we use nuclear matter obeying the APR 
equation of state [2lj and as a low density extension of the 
APR data we use p2l [23] . In order to apply our general 
results to other equations of state, away from chemical 
equilibrium we use the simple quadratic parameterization 
in terms of the symmetry energy employed in [24 . With 
the exception of ultra-heavy neutron stars the APR equa- 
tion of state allows only modified Urea processes |25fE7] 
which are the dominant microscopic processes that in- 
duce bulk viscosity. Direct Urea processes [271 [2S] are 
only possible in the core of a neutron star close to the 
mass limit which we study in addition. The dominant 
contribution to the shear viscosity in hadronic matter re- 
sults from non-Fermi liquid enhanced lepton-lepton scat- 
tering 

For the strange stars and the quark core of the hybrid 
stars we use quark matter obeying a simple equation of 
state in terms of parameters c, nig and B. 
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This is a generalization of the parameterization em- 
ployed in |30j from which the equilibrium pressure is 
obtained via the conditions of weak equilibrium /is = 
fJ-d = fJ'u+fJ'e and electrical neutrality. The shear viscosity 
in quark matter arises from non-Fermi-liquid-enhanced 
quark-quark scattering [31] and the bulk viscosity from 
non-leptonic flavor-changing weak processes |14l [52] . For 
our hybrid stars we make the assumption of local electri- 
cal neutrality, excluding the possibility of a mixed phase 
and its wealth of geometric structures. 

The equilibrium star configuration is then determined 
by the general-relativistic Tolman-Oppenheimer-Volkov 
(TOV) equations [33 . Characteristics of the examples 
that we consider are given in table |lj The susceptibilities 
B and C of the considered forms of matter are given in 
table [nj Finally, the reequilibration parameters f , S and 
X in the parameterization of the weak rate eq. ([T|) and 
the parameters of the shear viscosity eq. ([9| are given for 
the considered forms of dense matter in table IIIII Note 
that there are higher non-linearities in the weak rate of 
hadronic compared to quark matter, which in accordance 
with eq. Q yields a steeper rise of the viscosity with am- 
plitude so that suprathermal effects are even more im- 
portant in hadronic matter than in quark matter. 
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M[Mq\ 


M^ore [Mq\ 


R[km\ 


Uc [no] 


(n) [no] 


Q.K [kHz] 


NS 


1.4 


(1.39) 


11.5 


3.43 


1.58 


6.02 




2.0 


(1.99) 


11.0 


4.91 


2.46 


7.68 




2.21 


0.85 


10.0 


7.17 


3.37 


9.31 


SS 


1.4 




11.3 


2.62 


1.91 


6.17 




2.0 




11.6 


4.95 


2.43 


7.09 


HS 


1.4 


1.06 


12.7 


2.32 


1.17 


5.16 




2.0 


1.81 


12.2 


4.89 


1.84 


6.62 



Table I: Properties of the considered models of neutron stars 
(NS), strange stars (SS) and hybrid stars (HS). We show the 
mass of the star M, the mass of the core Mcore, the radius R, 
the baryon density at the center of the star Uc given in units 
of nuclear saturation density no , the average density (n) and 
the Kepler frequency Q.k- The neutron stars were obtained by 
solving the relativistic TOV equations for catalyzed neutron 
matter using the APR equation of state [21! with low density 
extension |22l I23| and the strange stars with a quark gas bag 
model with c = 0, iris = 150 MeV and a bag parameter B = 
(138 MeV)^. Large mass hybrid stars are only found when 
strong interaction corrections are considered, cf. |30| . and we 
find a 2 Mq star for c = 0.4, rUs = 140 MeV, B = (137MeV)'*. 



C. R-mode profile 



R-modes are normal modes of rotating stars and are 
obtained from a linear perturbation analysis around the 
static star configurations discussed above. We consider 
a star rotating with angular frequency f2. To simplify 
this analysis it is performed in a Newtonian approxima- 
tion and in a slow rotation expansion in (ttGnp), 
where p is the average density of the unperturbed star. 
The analysis of the damping due to bulk viscosity strictly 
requires an expansion to next to leading order. The os- 
cillation frequency depends on the frame. For an r-mode 
with angular dependence Y^_^_^ the frequency in a frame 
rotating with the unperturbed star ujr, which is relevant 
for microscopic quantities like the viscosity, and the cor- 
responding frequency Wi, observed by an observer in an 
inertial frame, are given by 
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where the ellipsis denotes further contributions involving 
next to leading order corrections of the potentials that 
are explicitly given in [2J. In this expression A denotes 
the inverse squared speed of sound 



A 



dp 

dp 



(13) 



evaluated at equilibrium. There are different conventions 
for the amplitude a in the literature and we follow the 
convention for the amplitude given in |34| but take into 
account the corrections to the latter result in [2]. This 
convention is usually used in the literature |48| and in 
this case the above expression breaks down for a>0 (1). 
For more details on the r-mode expression see Appendix 
[A} In general one must numerically solve a differential 
equation to obtain 6^q, which is the leading order cor- 
rection to the gravitational potential However, in 
the special case of a star with a constant density profile 
6^0 can be shown to be subleading compared to the first 
term in the inner parentheses of eq. (12l. For our nu- 



merical analysis below we will make the approximation 
to neglect the additional second order corrections in the 
slow rotation expansion given by the ellipsis in eq. (12 1, 



which amounts to replacing the Lagrangian perturbation 
by the Eulerian perturbation, but we include the second 
order corrections to the frequency eq. (111. General an- 



alytic results showed that this is a good approximation 
for the computation of the small amplitude instability re- 
gions Here we will give semi-analytic results for the 
saturation amplitudes, valid beyond leading order, that 
show the infiuence of the second order terms. 

We will see below that, in obtaining a precise assess- 
ment of the damping time of large amplitude r-modes, 
the radial dependence of the density perturbation plays 
a vital role. The radius enters eq. ( 12 ) explicitly and also 
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via the density dependence of the inverse squared speed 
of sound A (fig.[lj and the radial density dependence of 
the star (fig. 2 in [4 ). The radial variation of the den- 
sity is moderate in a strange star, but much more pro- 
nounced in neutron stars where the r-mode amplitude 
grows strongly in the outer regions of the star. 



where the function k has an analog expansion and reads 
to lowest order ko — 2/ (m+l). For dissipation via bulk 
viscosity the relevant quantity is the density compression 
due to the r-mode with dimensionless amplitude a. It 
vanishes to leading order in the slow rotation expansion 
and reads at next to leading order [5J [31] 



III. DAMPING TIME SCALES 

A. General expressions 

The amplitude of the r-mode oscillations evolves with 
time dependence exjp{iujt — t/r). We can decompose the 
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hadronic matter 



hadronic gas 
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A 



B 



n ' (4(1--23;)S)'^ 



(371 



N2/3 



3+ 



(1-c)m2 3(1-c)m2 



12(1- 



c 



4(l-2.)(n|f-f) 



(SjT^n 
6m , 



^273- 



3(1-c)m, 



Table II: Strong interaction parameters, defined in eqs. ([2| and (13 1, describing the response of the particular form of matter. 
In the case of interacting hadronic matter a quadratic ansatz in the proton fraction x parameterized by the symmetry energy 
S is employed. The expressions for a hadron and quark gas are given to leading order in n/m% respectively next to leading 
order in ms/jJ,. 



Weak process 



quark non-lcptonic 



hadronic direct Urea 
hadronic modified Urea 



f fMeV^^"'')] 



6.59x10-^2 



4.68x10"^^ (^^^ 



173 

T73- 



Xi 



1 

47r^ 



189 



3677r2 



X2 



21 



3677r4 



X3 



1835ir6 



Strong/EM process 


fj [MeV(^+")] 


a 


quark scattering 


1.98xl0V"^ (W^)^ 


5 

3 


leptonic scattering 


1.40x10^2 ( 


14 


5 
3 


nn-scattering 


5.46x10'' f " 


2 



Table III: Left panel: Parameters of the general parameterization of the weak rate eq. ([T]) for different processes of particular 
forms of matter which determine the damping due to bulk viscosity. The coefficients Xi parameterize the non-linear dependence 
on the chemical potential fluctuation /ia arising in the suprathermal regime of the viscosity which is relevant for large amplitude 
r-modes. Right panel: Parameters arising in the parameterization eq. Q of the shear viscosity for different strong and 
electromagnetic interaction processes. The leptonic and quark scatterings arise from a non- Fermi liquid enhancement due to 
unscreened magnetic interactions. 



decay rate 1/t as 
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Where tq, tb and ts are gravitational radiation, bulk 
viscosity and shear viscosity time scales, respectively. 
The time scale of the r-mode growth due to gravitational 
wave emission is given by [34] 
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and the damping time due to the shear viscosity can be 
written as ^ 



1 _ ^ (m - 1) (2m + 1) sii^A^+c^^R 



with the radial integral constants 
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where p is the energy density and Aqcd a generic QCD 
scale introduced to make the constant dimensionless. 
When several shells with inner radii i?,-'''and outer radii 



R^o^ and distinct phases and/or transport coefficients are 
present in a compact star, the damping time integral over 
the star decomposes into a sum of contributions of the 
individual shells The bulk viscosity damping time is 
given by [2J 
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Using the general expression for the bulk viscosity eqs. ^ 
and ^ and expressing the fluctuation in the conserved 
energy density Ap/p by the same fluctuation of the 
baryon density An/n, gives for the damping time the 
general expression 
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in terms of integrals over the individual shells. Defining 
the reduced density oscillation 
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Therefore we will in our numeric results given below ne- 
glect the contribution from the crust and only include 
the contribution from the core. The core does not have 
a sharply-defined boundary: we chose it conventionally 
to be at baryon density uq/A corresponding to the lowest 
point in the APR table, but check the dependence on this 
choice. 
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Figure 1: The density dependence of the inverse squared 
speed of sound A = dp/dp (which enters the r-mode profile 
multiplicatively) for the different forms of matter in table[ll|as 
well as generic polytropic models. The solid line represents 
interacting APR matter, the dashed line a hadron gas and 
the dotted line shows the result for a quark gas. The struc- 
ture at intermediate densities in the APR curve arises from 
phase transitions and the use of finite differences to compute 
the derivative, but due to the mild contribution of the denser 
inner regions of the star to the damping these, as well as 
the known problem that the APR equation of state becomes 
acausal at high density, have no influence on our results below. 



depending on only two independent parameters 

"~ Ka;i/<5 ^ (to + 1)k(<5+i)/« 



Note that in eq. ( 22 1 all local quantities can have differ- 



ent functional forms in different shells, as given in tables 
|TT] and |III| but to make the expression readable we do 
not show the explicit suffixes (1). Whereas strange stars 
are basically homogeneous and consist of a single phase, 
the crust of neutron and hybrid stars is extremely inho- 
mogeneous and complicated. Although there are no free 
protons and thereby no Urea processes, the ultra-heavy 
nuclei present in the inner crust as well as the clusters 
in potential pasta phases still feature analogous beta- 
processes. Since oscillations likewise push the system out 
of beta-equilibrium an analogous suprathermal enhance- 
ment of the bulk viscosity contribution from these phases 
is expected. However, there are to our knowledge no re- 
sults for the bulk viscosity in the inner crust, yet |35| . 



B. Approximate limits of the bulk viscosity 
damping time 

In the subthermal regime fiA ^ T the bulk viscosity 
eq. (|6| is independent of the r-mode amplitude, so the 
angular integral in eq. (19) is trivial. The damping time 



in the subthermal regime is then given by 
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in terms of the one dimensional radial integral 
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This expression was used to study the small amplitude 
instability regions in [4 . Here we want to study the large- 
amplitude saturation and therefore it is useful to obtain 
an analytic expression that includes the large-amplitude 
enhancement of the bulk viscosity. In the intermediate, 
linear regime and for / ^ 1 the general analytic approx- 
imation for the bulk viscosity eq. ([t]) is valid. Since this 
local condition has to be fulfilled everywhere in the star, 
the global parameters a and h must be smaller than cer- 
tain bounds that are determined by the particular prop- 
erties of the considered stars. We recall from ||T2| that the 
approximation is particularly useful for hadronic matter 
with modified Urea processes where it covers almost the 
entire range of physical local density amplitudes at mil- 
lisecond frequencies. A plot of the regions of validity of 
the individual analytic approximations of the bulk vis- 
cosity for different forms of matter is given in [19 . Anal- 
ogous to the low temperature/high frequency approxima- 
tion in the subthermal regime eq. ( 23 I , in the intermedi- 
ate, linear regime an explicit evaluation is possible. With 
the analytic expression for the bulk viscosity eq. (|7| the 
angular integrals over the spherical harmonics arising in 
the r-mode profile take the form 
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and this yields a result that apart from the evaluation 
of the remaining radial integrals is analytic (see also |13l 
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The dependence on all local parameters, like the equation 
of state, the weak rate, the density dependence of the 
particular star and its r-mode profile is contained in a few 
dimensionless radial integral constants {Aew is a typical 
electroweak scale) 
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but the dependence on the parameters of the r-mode evo- 
lution 51, a and T is entirely explicit in eq. (24l. The 
J = term in eq. ( 24 1 is precisely the approximate sub- 



thermal result eq. (23 1 in the considered approximation. 



The constants Vj = V2.j for the fundamental r-mode are 
given for several stars in table |IV[ Although these pa- 
rameters can vary significantly for different stars, it is 
quite striking that as far as the bulk viscosity is con- 
cerned, the complex details of the individual stars are 
encoded in a few constants. Note in particular that the 



parametric form eq. ( 24 1 remains valid for the full second 



order r-mode expression and only the constants eq. (25) 



are changed. At sufficiently large amplitudes the largest 
power in the sum in eq. ( 24 1 dominates and due to the 
connection 6 = 2N the bulk viscosity damping time be- 
comes temperature independent in this approximation. 
Note also that the integrals Vm^N, as well as the gen- 
eral expression eq. ( p2| feature an extremely pronounced 
radial dependence, both due to the explicit radial depen- 
dence and the radial dependence of the inverse squared 
speed of sound A, that strongly weights the outer parts 
of the star. 



C. Results for the damping times 

Using the expressions for the microscopic parameters 
given in tables a nd |III| in the general expressions eqs. 
(151, (16 1 and (20 1 we obtain the gravitational and vis- 



cosity time scales as a function of their dependent macro- 
scopic parameters. These are given in fig. [2] for the cases 
of a neutron star with damping due to modified Urea re- 
actions (left panel) and a strange star (right panel) by the 
solid lines as a function of temperature and for different 
amplitudes ranging from top to bottom from the sub- 
thermal result at infinitesimal amplitude to the extreme 
case a = 10. 

At sufficiently low temperature the strong increase 
of the bulk viscosity with the (local) amplitude An/n 



damps r-modes with large (global) dimensionless ampli- 
tude a at significantly shorter time scales. As found be- 
fore from the analytic expression eq. ( 24 ) , given by the 
dotted curves, the damping time is temperature inde- 
pendent in this low temperature and intermediate ampli- 
tude regime. In contrast, due to the generic form of the 
bulk viscosity, featuring a universal maximum, the damp- 
ing of large amplitude r-modes is not enhanced at high 
temperatures. As a direct consequence of the subther- 
mal maximum of the bulk viscosity [12^ , the correspond- 
ing "resonant" temperature where the damping time is 
minimal is at roughly 10^ K for strange stars and 10^^ 
K for neutron stars. Correspondingly, r-modes are en- 
tirely unstable at high temperatures. However, due to 
the strong suprathermal enhancement at low tempera- 
tures the damping undercuts the gravitational time scale 
at sufficiently large amplitude so that the r-mode growth 
will slow down and eventually saturate. The correspond- 
ing amplitudes are strikingly very similar for the two dif- 
ferent classes of stars, as will be discussed in more detail 
below. At very large amplitudes a^O (10) the damping 
times decrease again at all temperatures as a consequence 
of the behavior of the bulk viscosity [T^ • 

The dot-dashed curves in the neutron star plot on the 
left panel of fig. [2] show the damping time if the crust 
is assumed to start already at the higher density no/2 
instead of no/4 so that only the contribution from the 
correspondingly smaller core is taken into account [49]. 
Although the damping times are larger, as expected, the 
amplitude at which the viscous damping can saturate the 
mode is not drastically changed, so that our results given 
below remain qualitatively unchanged in this case. Ac- 
tually, when the damping from the crust would be prop- 
erly taken into account this should rather enhance the 
damping and decrease the r-mode amplitude compared 
to those obtained in this work. In the case of the strange 
star on the right panel of fig. [2] the dashed curves also 
show the analytic approximation discussed in appendix 
[B] where the star is assumed to be homogeneous. As can 
be seen, the corresponding expression eq. (B5) gives an 



approximate estimate for the damping time at all tem- 
peratures and for amplitudes up to the maximum of the 
viscosity, and only fails at higher amplitudes, where the 
bulk viscosity cannot saturate the r-mode anymore and 
where it is thereby not physically relevant. The devia- 
tions compared to the numeric result stem mainly from 
the fact that the density in the strange star is not entirely 
constant (fig. 2 in |3]). 

Previous neutron star analyses |13| l36] have employed 
an r-mode profile that does not feature the stronger addi- 
tional radial dependence due to the low density enhance- 
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Table IV: Radial integral parameters and static s atu ration amplitude of a m = 2 r-mode for the stars considered in this work. 
The constant J, S and Vi are given by eqs. (171, (18 1 and (251, respectively, using the generic normalization scales Aqcd ~ 1 
GeV and Aew ~ 100 GeV. Note that the subthermal parameter Vo corresponds to y in j4j where the subscript was omitted 
for simplicity and the strange star expressions are given to leading order in ms/jJ.. 



ment of the inverse squared speed of sound A, shown 
in fig.[lj for a realistic equation of state. E.g. the r- 
mode profile given in eq. (6.6) of f37] features roughly a 
generic r'^-dependence. The strong r-dependence in our 
present treatment, however, strongly amplifies the damp- 
ing in the outer regions of the star. Therefore, the con- 
tribution of the crust to the viscous damping should be 
relevant and would further decrease the saturation am- 
plitude. The current restriction of our analysis to the 
core presents therefore an upper bound for the satura- 
tion amplitude obtained when the damping of the whole 
star is considered. The second order effects in contrast 
increase the small amplitude instability region |2 and 
can correspondingly be expected to likewise increase the 
saturation amplitude. A more thorough treatment of all 
these effects in the future is clearly desirable. 



IV. SATURATION AMPLITUDES 

Because of the strong decrease of the viscous damping 
time due to the suprathermal enhancement the damp- 
ing can dominate at sufficiently large amplitudes. In this 
case the definition of the instability regions have to be 
extended. The latter are standardly defined in the sub- 
thermal regime and are independent of the amplitude. 
One could extend this concept by the definition of am- 
plitude dependent instability regions which would shrink 
with increasing amplitude. However, since the ampli- 
tude can neither be inferred from observation nor is it a 
parameter that can be dialed, but is rather determined 
dynamically, we refrain from this possibility and rather 
introduce the concept of a static saturation amplitude. 
The latter is defined by the amplitude at which the r- 
mode would saturate at fixed temperature and frequency 
and is given by the solution of the equation 



.(0 



(T) 4'^(r,r!,a,,t), 



(26) 



where I runs over the contributions from the different 
shells of the star. The boundaries of the above mentioned 
amplitude dependent instability regions are by defini- 
tion simply the contour lines asat{T, r2)=const. and the 
boundary of the classic instability region, in particular, 
corresponds to asat = 0. In case several solutions of 
eq. (26 1 exist, only the smallest one is physical and if no 



solution exists then viscous damping alone cannot sat- 
urate the r-mode according to this definition. Actually, 
at the same time the r-mode grows, the star generally 
also cools or reheats and spins down respectively up so 
that the above amplitudes do not have to be reached. 
In particular the star could leave the parameter regions 
where a saturation according to the above criterion is not 
possible before the r-mode can actually explode. 



A. Analytic approximation 

Similar to the analytic expression for the boundary of 
the instability region |3] an analytic expression for the 
static saturation amplitude can be obtained. The an- 
alytic linear approximation applies as long as the r- 
mode amplitude is small enough so that the bulk viscosity 
is sufficiently below its maximum. When the r-mode am- 
plitude is at the same time large enough that the highest 
power in eq. (24) dominates, the damping time becomes 



temperature independent. Both conditions are met suffi- 
ciently far inside the instability region so that the damp- 
ing time simplifies to 
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Figure 2: The relevant r-mode time scales for IAMq stars rotating at their Kepler frequency. Left panel: Neutron star. 
Right panel: Strange star. The dotted horizontal line presents the time scale tq associated to the growth of the mode due to 
gravitational wave emission. The dashed rising curve shows the damping time ts due to shear viscosity. The damping time 
tb due to bulk viscosity is given for different dimensionless r-mode amplitudes a = 0, 0.01, 0.1, 1 and 10 by the solid curves. 
The thin dotted curves correspond to the analytic linear approximation eq. (241 and are below the shown plot range for the 



largest amplitude. The thin dot-dashed curves on the left panel show the change when only a smaller core (ranging to a density 
of no/2 instead of no/4) is taken into account. The thin dashed curves on the right panel represent the approximate analytic 
expression eq. ( B5 1 given in the appendix which is not valid above the maximum of the bulk viscosity and therefore not shown 



for the large amplitude results. 
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At saturation this has to match the gravitational time scale 1 /tb + I/tq = which yields the general result 
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where k is defined by eq. (111. In the cases of strange 
stars with non-leptonic processes 5 = 2 and hadronic 
matter with modified Urea processes (5 = 6 this gives for 
the m = 2 r-mode 
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where J= J2 and Vi = ¥2.1 are given for the normaUza- 
tion scales used in table IV Here miso, A*300, Mia, Riq 
and ilms are the effective strange quark mass in units of 
150 MeV, the quark chemical potential in units of 300 
MeV, the stars mass in units of IAMq, the radius in 
units of 10 km and the angular velocity in units of 2tt 
kHz corresponding to a millisecond pulsar, respectively. 
For strange stars the above expression for the intermedi- 
ate amplitude bulk viscosity damping time has the same 
frequency dependence as the gravitational time scale eq. 
(151, so asat is basically constant throughout the insta- 



bility region. However, it rises with decreasing frequency 
for neutron stars where the frequency dependence of the 
bulk viscosity is weaker. In contrast to the analytic ex- 
pressions for the extrema of the instability region given 
in I?], which are very insensitive to the microscopic trans- 
port parameters, the saturation amplitude is more sen- 
sitive to the suprathermal bulk viscosity parameter V. 
Whereas the saturation amplitude of neutron stars still 
depends on V rather mildly due to the power 1/6, the 
saturation amplitude for strange stars obtained from the 
generic equation of state eq. ( 10 ) decreases with the "in- 



teraction parameter" c and even more strongly with the 
effective strange quark mass m^. 



B. Numeric solution 

Let us now discuss the numerical solution for the sat- 
uration amplitude. In the following plots figs.|3]to[7]the 
static saturation amplitude is shown as a function of tem- 
perature and amplitude and they feature generally 3 dis- 
tinct regions. The light (blue) surface shows the satura- 
tion amplitude where the r-mode growth is stopped by 
suprathermal damping. Due to the characteristic behav- 
ior of the bulk viscosity [12] which does not feature an 
amplitude enhancement for temperatures above the tem- 
perature Tmax eq. Q, the r-mode is not damped at all 
by viscous effects in the high temperature regime as de- 
noted by the dark (red) area on the right hand side. In 
the flat (green) region surrounding the instability region 
the r-mode is entirely stable and already damped by the 
shear or the subthermal bulk viscosity so that asat — 0. 

The left panel of flg.[3]shows the static saturation am- 
plitude for the m = 2 r-mode of a 1.4 Mq neutron star. 
Although this might be hard to see in certain regions of 
the plot, the saturation amplitude rises steeply, within 
a narrow interval, from zero at the boundary of the in- 
stability region towards its interior. At large frequencies 
it reaches a plateau value that is nearly independent of 
the temperature as predicted by the analytic expression 



eq. (27 1. As described by the latter expression, inside the 



instability region the saturation amplitude rises strongly 
with decreasing frequency and since the instability re- 
gion shrinks in width and eventually ends at low fre- 
quencies where the amplitude vanishes and the mode is 
damped, it features a peak- like structure. The maximum 
static saturation amplitude reached at the peak is in this 



case unphysically large whereas the plateau value at the 
Kepler frequency is still roughly 3.5. The suprathermal 
bulk viscosity can therefore in principle saturate r-modes 
within the lower part of the instability region. However, 
in the present case, where only the damping of the core is 
taken into account, the static saturation amplitudes are 
at the limit where a standard r-mode analysis is valid. 
Moreover these amplitudes are so far larger than those 
of alternative saturation mechanisms f71-fTCT|. It is inter- 
esting to mention once more, though, that the extreme 
radial dependence of the r-mode profile eq. ( [T2| strongly 
weights the outer regions of the star due to power law 
dependences of the inverse bulk viscosity damping time 
eq. ( 19 1 with exponents O (30) for neutron and hybrid 
stars which is further enhanced by the density depen- 
dence of the inverse speed of sound. The contribution 
of the crust could thereby be decisive to obtain a realis- 
tic estimate of the impact of the non-linear viscosity. In 
this context it is also important that a similar enhance- 
ment of the bulk viscosity has been found for superfluid 
matter [3S] . As noted in there is a second instability 
region at high temperatures above 10^^ K and as argued 
above the suprathermal bulk viscosity cannot saturate 
the r-mode in this high temperature regime. It is an in- 
teresting question if the r-mode can become large during 
this initial part of a star's evolution, and if so whether 
the r-mode is saturated at sufficiently small values by 
other non-linear mechanism or if the r-mode growth is 
not stopped before it reaches the regime where the struc- 
tural stability of the star is at stake. In the latter case 
this instability phase might extend the violent supernova 
stage and actively shape the remnant by additional mass 
shedding and thereby determine its initial size and angu- 
lar momentum. We will discuss these points further in 
the conclusion. 

The saturation amplitude for the IAMq strange star 
is shown on the left panel of fig.|4] Since the maximum of 
the stability window is above the Kepler frequency there 
are in this plot two separate parts of the instability re- 
gion. As predicted by eq. ( p7| the plateau value of the 
saturation amplitude in the lower part is approximately 
temperature and frequency independent. Strikingly it is 
of similar size as the saturation amplitude for the IAMq 
neutron star at its Kepler frequency. The high temper- 
ature part of the instability region where the viscosity 
again cannot saturate the r-mode extends in this case 
down to lower temperatures than for neutron stars. 

The left panel of fig. |5] shows the corresponding plot 
for the 1.4 Mq hybrid star. As found previously in the 
instability region has here three parts that are separated 
by two stability windows arising from the resonant be- 
havior of the bulk viscosities in the different shells. Due 
to their parametrically different temperature dependence 
the bulk viscosity of the quark shell dominates at low 
temperature, whereas the bulk viscosity of the hadronic 
shell dominates at high temperatures. Correspondingly 
the saturation amplitude in the low temperature part 
of the instability region shows the qualitative behavior 
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Figure 3: The static saturation amplitude, at which the r-mode growth is stopped by suprathermal viscous damping for the 
APR neutron stars. Left panel: IAMq . Right panel: 2.0 Mq. The light (green) shaded area denotes the stable region where 
the r-mode is damped away. At large frequencies a plateau with amplitudes 0(1) is reached. In the dark (red) region at high 
temperatures the r-mode is entirely unstable and cannot be saturated by viscous effects. 




Figure 4: The saturation amplitude for the considered strange stars. Left panel: IAA^q . Right panel: 2.0 Mq . In the latter 
case the suprathermal viscosity cannot stop the r-mode instability at frequencies larger than the majdmum frequency of the 
stability window (where the saturation amplitude diverges) - for the considered star slightly below the Kepler frequency - as 
well as in the high temperature part of the instability region. The saturation amplitudes of the plateau in the lower part of the 
instability region are of the same order as in the hadronic case shown in fig.[3] 



found for strange matter whereas the intermediate tem- 
perature part shows the quahtative behavior found for 
hadronic matter. Since the region where the peak in 
fig. [3] is located is "cut out" by the stabiUty window, the 
remaining peak of the hadronic intermediate part of the 
instabihty region in fig. [5] reaches only a much lower am- 
plitude. 

The result for the heavy 2.0 Mq neutron star is given 
on the right panel of fig.|3] As had been found previously 
in lU the instability region is larger for such heavy stars. 
The figure shows in addition that saturation occurs at a 



somewhat higher amplitude. The 2.0 Mq strange star 
is given on the right panel of fig.|4j In this case the 
maximum of the stability window is below the Kepler 
frequency. Similar to the high temperature behavior dis- 
cussed before, the r-mode cannot be damped by viscous 
effects above this maximum. It is interesting to recall 
from jl] that in the case of quark matter an approximate 
analytic expression for the location of the maximum of 
the stability window exists 
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Figure 5: The saturation amplitude for the considered hybrid stars. Left panel: 1.4 . Right panel: 2.0 Mq . The saturation 
in the low temperature part of the instability region is mostly established by the bulk viscosity of the quark core, whereas the 
saturation in the mid temperature part comes mainly from the hadronic shell. 
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where f = T / fj,^. This shows that m addition to a large 
star mass, a small effective strange quark mass in the 
quark matter equation of state eq. (10 1 increases the to- 



tal instability region both at high frequency and high 
temperature. In contrast, for the heavy 2.0 hybrid 
star shown on the right panel of fig.jS] such a total in- 
stability region does not arise since although the quark 
core cannot saturate the r-mode, the hadronic shell alone 
still provides sufficient damping to do so. In summary r- 
modes in massive stars are more unstable than in light 
stars since both their small amplitude instability regions 
are larger and they are less efficiently saturated by the 
large amplitude enhancement of the bulk viscosity. 

The left panel of fig. [6] shows the static saturation am- 
plitude for a neutron star with 2.21 Mq which is the max- 
imum mass allowed by the APR equation of state. In this 
case direct Urea reactions are possible in a small inner 
core region of mass 0.85 M©. As had already been ob- 
served in [31 , direct Urea reactions only slightly alter the 
instability boundary by a small notch at its right hand 
side. Since suprathermal damping from outer layers dom- 
inates due to the strong radial dependence of the r-mode 
the static saturation amplitude is likewise only slightly 
reduced by the small direct Urea core. However, because 
the size of the inner direct Urea core depends strongly 
on the equation of state and there are equations of state 
where the direct Urea core is considerably larger, we show 
on the right panel of fig.|6] for comparison the (unphys- 



ical) case that the direct Urea reactions are artificially 
switched on in the entire core. This represents an up- 
per limit for the possible effect of direct Urea reactions 
and shows that in this extreme case the static saturation 
amplitude at large frequency is reduced and the increase 
towards lower frequencies is considerably weakened ac- 
cording to the 1 /ri-behavior predicted by eq. ( 27 1 . 

In contrast to the previous results that evaluated the 
damping time eq. ( |19| numerically the top, left panel 
of figJT] employs the approximate analytic expression 
eq. (|24[) for the m = 2 mode of the IAMq neutron star. 
Comparing it to the numerical result in fig. |3] shows that 
the corrections are very small and because the maxi- 
mum of the bulk viscosity of hadronic matter with mod- 
ified Urea reactions is reached only for large amplitudes, 
eq. (24 1 provides a very good approximation in this case. 



In contrast, the use of the linear approximation which 
neglects the large amplitude decrease of the bulk viscos- 
ity, strongly overestimates the damping for the case of 
strange stars and misses the previously discussed total 
instability region at high frequency in fig.^ 

Fig. [7] also shows the saturation amplitudes of differ- 
ent multipole r- modes, given for the first four multipoles 
TO = 2 to 5 of the 1.4 Mq neutron star . The higher 
multipoles saturate at lower amplitudes than the to = 2 
and therefore the use of the linear approximate is well 
justified in this case. Interestingly, although the right 
segments of the lower part of the instability boundary of 
these higher order r-modes had recently been shown to 
be very similar to that of the fundamental to = 2 mode 
[4J, fig. [t] shows that although the peak value of the sat- 
uration amplitude of these modes decreases, the value 
at the Kepler frequency stays nearly constant. There- 
fore, these higher multipoles could be relevant for the 
spin-down evolution since the spin-down torque due to 
gravitational wave emission depends strongly on the am- 
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Figure 6: The static saturation amplitude, at which the r-mode growth is stopped by suprathermal viscous damping for the 
APR neutron stars at the maximum mass 2.21 Mq, where direct Urea processes become allowed. Left panel: direct Urea is 
only allowed in a small inner core region, see Table [I] Right panel: the same model when direct Urea is artificially turned on in 
the entire core. 




Figure 7: Saturation amplitudes for the first four multipole r-modes of the 1.4 Mq neutron star (top, left: m — 2; top, right: 
m = 3; bottom, left: m = 4; bottom right: m = 5). The results are obtained in the linear approximation eq. (24 1. 
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plitude [S]; for sufficiently small amplitude modes this 
dependence is quadratic. So if the suprathermal damping 
is responsible for the r-mode saturation, the restriction 
to the lowest order mode, that had been employed in all 
present analyses, should present only a first approxima- 
tion. 

The saturation amplitudes of the different 1.4 Mq stars 
are finally compared with each other and the analytic 
approximation eq. (24l in fig.js] Surprisingly, all stars 
feature saturation amplitudes of the same order of mag- 
nitude for millisecond pulsars, despite their very differ- 
ent microscopic and structural aspects. As noted before, 
for larger oscillation periods hadronic stars and to some 
extent also hybrid stars feature considerably larger satu- 
ration amplitudes than strange stars due to the paramet- 
rically different frequency dependence, see eq. (27 1. The 



analytic approximation yields in most cases a reasonable 
approximation to the full results with errors below the 
10% level. In general the analytic result overestimates 
the actual amplitude since it only describes the result far 
away from the boundaries and boundary effects play a 
role. In contrast at large frequencies the analytic approx- 
imation underestimates the saturation amplitude since 
the considered frequencies are already close to the criti- 
cal values, where the saturation amplitude diverges, see 
fig. Q. Nevertheless, the analytic approximation pro- 
vides an important and reliable estimate for the order of 
magnitude of the static saturation amplitude which, as 
will be discussed in more detail below, provides an upper 
limit for saturation amplitudes taken in dynamical star 
evolutions. 



V. CONCLUSIONS 

Using the recent general results for the bulk viscosity 
that include its non-linear behavior at large amplitudes 
we have derived expressions for the r-mode damping time 
that show that in the regime below the resonant tem- 
perature of the bulk viscosity, large amplitude r-modes 
are damped on considerably shorter time scales than low 
amplitude oscillations. In contrast, the universal maxi- 
mum of the bulk viscosity found in [T^] implies that at 
very high temperatures and frequencies r-modes cannot 
be damped at all by viscous effects since there is no en- 
hancement in the suprathermal limit. We find that for 
most stars considered in this work the corresponding crit- 
ical frequency is above the Kepler frequency. On the 
other hand the r-modes of all considered stars are un- 
stable at temperatures that are expected to be present 
when a proto-neutron star is created. At lower temper- 
atures our results lead to an extension of the concept of 
the instability region of an r-mode since the latter is only 
initially unstable at small amplitudes but the suprather- 
mal viscous damping can saturate the r-mode growth at 
finite amplitudes. We find that well within the instabil- 
ity region the static saturation amplitude agat defined 
in the text is temperature independent and takes values 



0(1) at milli-second frequencies for all considered stars. 
This is incidentally the order of magnitude that had been 
assumed in early r-mode analyses |S]. Yet, the static val- 
ues obtained here represent only an upper limit for the 
actual amplitude reached in the dynamic evolution. Our 
numeric results are confirmed by approximate analytic 
expressions which reveal the dependence of these results 
on the various underlying parameters. We also studied 
higher multipoles and find that although the first few 
multipoles have instability regions that are sizable, they 
feature similar saturation amplitudes as the fundamental 
m = 2 mode for millisecond pulsars and could thereby be 
relevant. 

It is interesting to compare our saturation mechanism 
and the obtained results for the saturation amplitudes 
with previously proposed mechanisms. In general when 
there are different competing saturation mechanisms, the 
one with the smallest saturation amplitude should dom- 
inate and effectively saturate the mode. Explicit numer- 
ical analyses of the general relativistic hydrodynamical 
equations [51 [51 would present the ideal way to study 
the saturation and star evolution. Whereas some of these 
studies find saturation only at large amplitudes, in oth- 
ers the r-mode can be completely destroyed by the decay 
into daughter modes once it exceeds amplitudes O (10~^) 
|10| . see also [23 • However, the numerical complexity 
limits these analyses so far to unphysically large values 
of the radiation reaction force that are orders of magni- 
tude above the physical value and it is not clear to what 
extent the obtained results can be extrapolated to the 
physical case. Another proposed saturation mechanism 
relies on the non-linear coupling of different oscillation 
modes [7] |H1 [Ml ilQj . These analyses find that this mech- 
anism could saturate r-modes at amplitudes as low as 
O (lO~^). Due to the considerable difficulties of a com- 
plete description of such a mode coupling mechanism, 
these analyses have to rely on model systems of generic 
coupled oscillators without a detailed connection to the 
complicated coupling of collective star oscillations. In 
summary, within the present approximation to neglect 
the neutron star crust, competing mechanisms will very 
likely dominate and saturate the r-mode at lower val- 
ues than the suprathermal enhancement of the viscosity. 
However, these mechanisms still involve simplifications 
and uncertainties. Our novel saturation mechanism, in 
contrast, relies on standard viscous effects and micro- 
scopic physics that is quantitatively well understood. 

Let us now discuss the implications of our results for 
the spin-down of compact stars. In the supernova for- 
mation process where a much larger star contracts to a 
very compact object that takes over the angular momen- 
tum it seems plausible that fast rotating proto-neutron 
stars could be formed which spin with frequencies close 
to the Kepler limit. According to our results r-mode os- 
cillations are unstable in this initial hot stage T>10^'^K 
and cannot be saturated by viscous effects for all consid- 
ered forms of dense matter. Generically, the cooling is 
very fast in this regime so that the evolution could leave 
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Figure 8: Comparison of the saturation amplitudes for the different 1.4 Mq stars. Left panel: Stars spinning with a period of 
1 ms. Right panel: Same for stars rotating with a period of 4 ms. Shown are the considered neutron star (solid), the hybrid 
star (dashed) and the strange star (dotted). The thick curves present the numerical results and the thin horizontal segments 
denote the analytic values obtained from eq. (27 1. 



this instability region before large amplitude r-modes de- 
velop or spin down the star. The star will then cool until 
it reaches the lower instability zone and the r-mode devel- 
ops. According to fig.[8]in this regime the r-mode can be 
saturated by viscous damping. For strange stars such a 
saturation does not seem to be required at all since the in- 
stability region is in this case located at comparably low 
temperatures [41' where the cooling becomes slow and 
the star either quickly spins down [42^ . or when reheat- 
ing effects are considered it reheats again [JS], and leaves 
the instability region before the amplitude becomes large. 
In this case the evolution wiggles around the instability 
line thereby spinning down the star, but this can take 
billions of years due to the strong reheating. 

In contrast in the case of neutron and hybrid stars 
without strangeness, the instability region is reached at 
large temperatures where cooling is still fast and reheat- 
ing effects are moderate, so that the evolution quickly 
penetrates the instability region and a saturation mech- 
anism is indeed required to stop the r-mode growth [5J. 
Since the static saturation amplitude increases continu- 
ously at the boundary of the instability region the dis- 
cussed static value does not have to be reached but a 
dynamic equilibrium could be established at a lower sat- 
uration amplitude that is reached once the r-mode is suf- 
ficiently large that the spindown becomes efficient. Due 
to this the viscous saturation could dominate competing 
saturation mechanisms. Once the r-mode is saturated, 
the question is which one of two competing processes. 



cooling or spin-down, is faster. Since the cooling is slow- 
ing down at lower temperatures it is likely that the spin- 
down wins and the evolution leaves the instability region 
near its lower boundary. In this case no young compact 
stars with frequencies larger than a tenth of the Kepler 
frequency would be possible which is in good agreement 
with observations. An answer to the above questions re- 
quires a detailed study of the combined spin-down and 
cooling evolution of the star which will be presented else- 
where. 

Strikingly our results suggest even another possibil- 
ity for the spindown of young stars that would be even 
faster and more violent. The core bounce during the 
supernova process should excite rather large amplitude 
oscillation modes in the forming compact core. Since r- 
modes are unstable in this regime [3] these will grow fur- 
ther. Because of the initial high temperatures, neutrinos 
are trapped inside the proto-neutron stars for roughly a 
minute |44]. Since a neutron star crust, that could pro- 
vide an efficient damping mechanism [11' , is not formed 
at this point and as our results show viscous effects can- 
not stop the r-mode growth, the amplitude could indeed 
become large if other non-linear saturation mechanisms 
likewise cannot operate efficiently in this turbulent en- 
vironment. In this case the loss of angular momentum 
could proceed not by gravitational wave emission but by 
actual mass shedding and thereby effectively as an ex- 
tension of the supernova explosion that is driven by r- 
modes. Since such a violent spindown should be fast 
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the star could end up at the lower boundary of the high 
temperature instability region before the star becomes 
transparent to neutrinos and the cooling process starts. 
Clearly, in this initial stage, which cannot rigorously be 
separated from the aftermath of the supernova explosion, 
the dynamics is highly non-linear and our simple r-mode 
analysis might not directly apply. Whether such a mech- 
anism is feasible will therefore require further study, but 
this mechanism would naturally explain the observed ab- 
sence of fast, young pulsars independent of their internal 
composition and it is striking that the frequencies of the 
high temperature instability boundary also seem to agree 
well with fastest pulsars that are young enough that they 
cannot be spun up by accretion |45j . 

Finally, r-modes should also be relevant for old accret- 
ing stars in binary systems that are spun up and could 
enter the instability region at low temperatures from be- 
low [46l. As discussed in [4j, strange and hybrid stars 
feature stability windows at low temperatures where the 
r-mode is absent, so that such stars could accelerate to 
frequencies close to the Kepler frequency. In contrast for 
neutron stars there is no stability window at low temper- 
atures so that an accreting star would enter the unstable 
regime already at low frequencies. Recall that the satura- 
tion amplitude of neutron stars due to bulk viscosity has 
a characteristic form with a pronounced peak close to the 
minimum of the instability region. In case the r-mode is 
saturated by suprathermal bulk viscosity, the steep rise 
of the amplitude close to the maximum should spin down 
the star quickly and so that it cannot penetrate deep into 
the instability region. This means that such stars should 
cluster close to the boundary which might be a signa- 
ture once more observational data for the temperature of 
compact stars becomes available. 
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can be solved analytically by matrix inversion. The ex- 
pression for SU reads to leading order in a slow rotation 
expansion 



SU 



TT {m + if {2m + ly. 



/ r \ m+l 



(Al) 

The velocity fluctuation is obtained from 5U by applica- 
tion of a differential operator [47J and yields to leading 
order in a slow rotation expansion |34| 



r) 

This expression provides the definition of a within the 
convention used here. In spherical coordinates this ex- 
pression yields the explicit form 



2™m! 



An (m + 1) 



[r) 



m—l i(mip+LLj 



i9 + cos ( 



where and (j) are unit vectors in polar and azimuthal 
direction. With this definition and for a — 1 the max- 
imum value, taken at the equator and in direction of 
9, is roughly Sv/v w 0.3, so that the approximation 
breaks down for a 3> 1 since mass shedding will occur for 
fast spinning stars. The corresponding maximum density 



fiuctuation 6n/n obtained from eq. ( 12 1 is more than an 
order of magnitude smaller. 

In contrast in [2^ an alternative convention a' of the 
amplitude was introduced that is related to the above a 

by 



a 



' TT {m + ly {2m + 1)1 , 

a 

m 



i.e. defined by eq. (All without the algebraic prefactor. 



Appendix B: Approximate result for strange quark 
matter 



Appendix A: R-mode and amplitude conventions 

In this appendix we review a few standard expressions 
for r-modes and discuss different conventions used in the 
literature. The form of the r-mode oscillation is most con- 
veniently derived [47 from the underlying equations that 
determine the fluctuation of the potential 5U = (5/i + (5$, 
where h is the enthalpy and <i> the gravitational poten- 
tial, since then the hydrodynamic Euler equation for the 
harmonic, cylindrically symmetric perturbation reduces 
from a differential to an ordinary linear equation and 



Whereas the bulk viscosity of semi-leptonic processes 
in hadronic and quark matter require in general a nu- 
meric solution, for the dominant contribution from non- 
leptonic processes in strange quark matter an approxi- 
mate solution valid in both subthermal and suprathermal 
regimes is possible. The approximate analytic result for 
the bulk viscosity in the suprathermal regime obtained 
from a Fourier analysis, is given by [T5] 



2 ^/ 973xrBC^r'^-^ / Ara. ^ 
3^3 \ ^ uj V 



(Bl) 
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where x = Xi ^iid 



= ^ ^y^^-zj'+(yi2^+zj'-2j (B2) 

A very good parameterization valid in both the sub- and 
suprathermal regime is given by [12 



(B3) 



Since the functional form of the suprathermal viscosity 
eq. (B2) is still complicated and does not allow to per- 



form the necessary subsequent integrations to obtain the 
r-mode damping time in an analytic form, we give an 
approximate analytic result that is valid up to the max- 
imum of the bulk viscosity. To this end we perform a 
global polynomial interpolation to the function h{z) in 
the interval [0,Zmax]- In order to appropriately describe 
the low amplitude behavior and the qualitative form be- 



low the maximum requires at least a quartic polynomial 
which is then uniquely determined as 



hpoi{z) = z- 



2V3 



z' + 



27 



324v^' 



(B4) 



The leading linear term in eq. ( B4 ) reproduces the ap- 



proximate intermediate linear result given by Madsen 
[14], whereas the other terms ensure the proper large 
amplitude saturation. The analytic form and the poly- 
nomial approximation agree in the relevant region below 
the maximum point-wise on the 5% level and the corre- 
sponding integrals required for the damping time to even 
better accuracy. 

The density in a strange star is nearly constant and so 
the density dependent quantities can be approximated by 
their value at the radius of the star, denoted by the suffix 
R. Performing the integration over the r-mode profile 
eq. (12), we find the approximate analytic result for the 



viscous damping time 



(B5) 

I 

where g is the polynomial 

, , 151875 2 1063125 . 290631796875 4 

q(x) — X X H X X 

' 9044 7429 587723048 
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